arXiv:1508.03623v2 [cond-mat.mes-hall] 22 Oct 2015 
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We study a proximity junction between a conventional s-wave superconductor and a conductor 
with Rashba spin-orbit coupling, with a specific focus on the spin structure of the induced pairing 
amplitude. We find that spin-triplet pairing correlations are induced by spin-orbit coupling in both 
one- and two-dimensional systems due to the lifted spin degeneracy. Additionally, this induced 
triplet pairing has a component with an odd frequency dependence that is robust to disorder. 
Our predictions are based on the solutions of the exact Gor’kov equations and are beyond the 
quasiclassical approximation. 


I. INTRODUCTION 


The generation of triplet superconducting correlations 
via proximity to a conventional superconductor (S) is a 
topic that has attracted a lot of attention recently in 
the context of superconducting spintronics,lill^ whereby 
triplet Cooper pairs with spin projection ±1 play the 
same role as electron spins in conventional spintronics.*^ 
Most of the research in this area, both experimental and 
theoretical, has focused on proximity junctions involv¬ 
ing ferromagnets (F), which induce triplet correlations 
by lifting the spin degeneracy (for a comprehensive re¬ 
view, see Refs. [9] and [in]). The induced triplet correlations 
are also of fundamental interest because they have an 
odd frequency dependence; by the Pauli principle, these 
triplet states must be isotropic in momentum and there¬ 
fore robust to disorder 

Because spin-orbit coupling (SOC) also lifts the spin 
degeneracy, one might expect odd-frequency triplet pair¬ 
ing to emerge when a ferromagnet is replaced by a spin- 
orbit material as well. That the pairing symmetry of a 
bulk spin-orbit-coupled superconductor is known to be a 
mixture of singlet and tripleipE^ only enhances these 
expectations. However, Liu et al. (Ref. [T7|) showed, 
by solving the exact Bogoliubov-de Gennes (BdG) equa¬ 
tions, that the proximity-induced pairing amplitude in 
a ID metal with Rashba-typ^^^ SOC (R) has no triplet 
component. Previous studies of proximity junctions in¬ 
volving both SOC and ferromagnetism that were carried 
out to leading order in the quasiclassical approximation 
also did not find any triplet pairing induced by SOC.^^ 
Only very recently has the singlet-triplet mixing effect 
of SOC in such structures been note d by working be¬ 
yond leading quasiclassical order.l^^HI] However, because 
the induced triplet pairing only shows up to second or¬ 
der, this is expected to be a very weak effect in any 
materials for which the quasiclassical methods are ap¬ 
plicable. Spin-orbit scattering confined to the interface 
has been shown to generate triplet pairing in 3D bal¬ 
listic superconductor/normal-metal junctions; however, 
this triplet component is anisotropic in momentum and 
therefore sensitive to disorder.l^ To date, odd-frequency 
triplet pairing induced by SOC in the proximity effect 
has not been explicitly investigated. 


In this paper, we show by directly solving the fully 
quantum-mechanical Gor’kov equations that spin-triplet 
superconducting correlations are induced by Rashba SOC 
in both ID and 2D proximity junctions. However, we find 
that the induced triplet component in ID vanishes when 
integrated over the momentum; this result is in agree¬ 
ment with Ref. [HI In 2D, we show that the induced 
triplet amplitude has an odd-frequency component that 
is isotropic in momentum. In agreement with Ref. 1191 
we also find that the triplet pairing induced by SOC 
vanishes to leading order in the quasiclassical approxi¬ 
mation. For this reason, our results are most relevant to 
materials that have a spin-orbit energy scale comparable 
to the Fermi energy,!^ when no quasiclassical expansion 
can be made. Examples include the surface states of no¬ 
ble metal^^ and semi-metallic bismuth,!^ InSb quantum 
wires,!^ as well as the bullP^ and surfac^^ states of the 
bismuth tellurohalides. The proximity effect in materials 
with strong SOC is also relevant to recent experiments 
probing the existence of Major ana fermions in semicon¬ 
ductor quantum wires.^^ 


The remainder of the paper is organized as follows. In 
Sec. |T^ we demonstrate that the presence of triplet pair¬ 
ing and the spatial dependence of the induced pairing 
amplitude can be deduced by considering Andreev re¬ 
flection processes. We review the well-known S/F prox¬ 
imity effect from this point of view in Sec. ra before 
discussing the qualitative similarities and differences of 


the S/R proximity effect in Sec. HB In Sec. HI 


we 


solve the Gor’kov equations in the S/R geometry to show 
that triplet pairing is induced via the proximity effect. 
Our explicit model is described is Sec. HI Aj and de¬ 
tails of the calculation are given in Sec. jlllB Solution 


details for the superconductor and Rashba regions are 
given in Sec. EIIB ^ and Sec. |HI|B.2[ respectively, while 
the enforcement of boundary conditions is discussed in 
Sec. IHIIB.3I Results in ID and 2D are discussed in 
Secs. HIC and HID In Sec. [HIE we show that the in¬ 
duced triplet pairing amplitude vanishes to leading order 
in the quasiclassical limit. Our conclusions are given in 
Sec. El 
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FIG. 1. (Color online) Bogoliubov-de Gennes excitation spectra of (a) a ID ferromagnet and (b) a ID Rashba metal with 
Rashba vector g = (0,0, —ifcj;). Colors denote different spin states, while solid (dashed) lines represent electron (hole) bands. 
An incident electron (e^ or 64 .) with energy e < A is Andreev reflected as a hole of opposite spin {hi or h^-, respectively). The 
presence of a triplet pairing component with zero spin projection results from the lifted spin degeneracy, (c) Fermi surfaces 
of a 2D Rashba metal, with Rashba vector g = {iky, —ikx,Q), taking Ef > 0. Because the two Rashba subbands are not 
associated with a definite spin, an incident electron (ei or 62) with momentum \ky\ < kpi can be Andreev reflected into either 
subband {hi or ^, 2 ). The total pairing amplitude is a sum of oscillatory and nonoscillatory terms, resulting from interband and 
intraband Andreev reflection, respectively. Note: though they appear different, the ID picture displayed in (b) and the ky = 0 
limit of the 2D picture displayed in (c) are physically equivalent. 


II. TRIPLET PAIRING FROM ANDREEV 
REFLECTION 

The proximity effect can be understood qualitatively 
through Andreev reflection, whereby an electron (hole) 
with excitation energy e below the superconducting gap 
A is retroreflected as a hole (electron) with opposite 
spin, thus injecting a spin-singlet Cooper pair into the 
superconductor.!^ While typically studied in the context 
of transport,!^il®! Andreev reflection can also give insight 
into the type of pairing that is induced by the proximity 
effect. Due to the coherence between electrons and holes, 
it is these processes which allow the pairing amplitude to 
penetrate into the nonsuperconducting material. 

By examining the eigenstates of the BdG equation,!^ 
we will show that SOC on its own is sufficient to induce 
triplet pairing in the proximity effect. For the sake of 
completeness, we will first review the case of triplet pair¬ 
ing in S/F junctions before extending the discussion to 
S/R junctions. 


A. Triplet pairing in S/F junctions 

We hrst present a brief review of the proximity effect in 
ferromagnetic materials, following qualitative arguments 
similar to those given in Ref. 1361 In a ferromagnet, the 
BdG equation reads 

(’'“"o''''" + = (') 

where Hq = {—d^. E k"^^)l2m — Ep, with /cy being the 
conserved momentum along the S/F interface, J is the 


ferromagnetic exchange field, is a Pauli matrix (we 
denote 2x2 matrices in spin space by a hat), and e) 
is a spinor wave function in Nambu G spin space. The 
ID BdG excitation spectrum, containing spin-split bands 
with distinct Fermi momenta = \/‘2m{Ep ± J) Ki 

kp ± J/vp, is shown in Fig. [^a). 

Because an incident electron [denoted by or in 
Fig. I^a)] must be Andreev reflected as a hole of oppo¬ 
site spin {hi or h^, respectively), the resulting Gooper 
pairs in a ID ferromagnet acquire a finite center-of-mass 
momentum Q = kp^ — kpi « ^hjvp, thus causing the 
Cooper pairing amplitude to oscillate in space with a 
period TTVp/h. [For arbitrary incidence, as is possible 
in higher dimensions, the center-of-mass momentum is 
modified to Q = 2J/vp cos9, where 0 is the angle the 
momentum makes with the interface normal.] The oscil¬ 
lations of the pairing amplitude are a direct consequence 
of the broken time-reversal symmetry, which implies that 
E\f\ yf A_k^. As a result, paired states with opposite 
spins must have a finite total momentum. These oscilla¬ 
tions are responsible for several interesting phenomena, 
including the possibility of a tt phase shift of the current- 
phase relation in S/F/S Josephson junctions, as well as 
nonmonotonic variations, as a function of the thickness of 
the ferromagnetic layer, of both the density of states and 
the superconducting critical temperature of S/F multi¬ 
layers (see Refs. [S] and [TO] and the references therein for 
more detail). 

As is more relevant to the current work, the lifting of 
the spin degeneracy (in this case by the exchange field) 
also results in the generation of a triplet pairing com¬ 
ponent that has zero spin projection. This is true for 
any material in which spin is an eigenstate of the Hamil¬ 
tonian and can be seen particularly easily in ID junc¬ 
tions by considering the contributions of the > hi 












3 


and —)■ Andreev reflection processes to the pairing 
amplitude. In terms of the BdG wave functions Ua{x,e) 
(electrons) and Va{x,£) (holes), the triplet component of 
the pairing amplitude with zero spin projection is 


J't (x, X ^ s') 2 X , s) A (x, X , 5 ) 


v^{x, e)u1{x', e) + vi{x, e)ut{x', e). 


( 2 ) 


The first term results from the —>■ process, while 

the second term results from the —>■ process. While 
the most general case is treated in Appendix]^ we can 
express the triplet pairing amplitude in the case where 
the electrons are perfectly Andreev reflected by 


J'j{x,x',e) ~ a^(£)e**^T'*^e 




(3) 

where is the momentum of a right-moving (left- 

moving) spin-up (spin-down) electron (hole) and a^(^)(£) 
is the amplitude for Andreev reflection from a spin-down 
(spin-up) electron to a spin-up (spin-down) hole. As 
shown in Appendix]^ the Andreev reflection amplitudes 
are related by a^(£) = —a^(£). Therefore, 


J^l{x,x',£) 






(4) 

Because this quantity can only vanish when 

any lifting of the spin degeneracy will result in a nonzero 

triplet pairing amplitude. 


B. Triplet pairing in S/R junctions 


While the proximity effect in S/R junctions is very 
similar to that in S/F junctions, there are several impor¬ 
tant qualitative differences we will now discuss. In the 
presence of SOC, the BdG equation is given by 


^ Ho — iAg • a 0 

V 0 —iio — iXg ■ a* 


ip{x, e) = £tjj{x, e), 


(5) 

where A is the SOG constant and g is the spin-orbit vec¬ 
tor, which satisfies g(k) = —g(—k). 

In a strictly ID system, where there is no well-defined 
spin quantization axis, the two representations of the 
Rashba vector g = {0,—ikx,0) and g = {0,0,—ikx) 
are physically equivalent. To discuss Andreev reflec¬ 
tion in the ID case, we choose the latter representation, 
so that spin remains an eigenstate of the BdG Hamil¬ 
tonian. The BdG excitation spectrum, which contains 
two Rashba subbands with Fermi momenta kpi(^ 2 ) = 
(mA)^ + 2mEp A rn\, is shown in Fig. B>>) . While 
we take Ep > 0 in Fig. 0 b) , so that the Fermi en¬ 
ergy lies above the Dirac point, the following physical 
arguments are equally valid \i Ep < 0. Because each 
Rashba subband can be associated with a definite spin, 
an electron incident on the S/R interface must be An¬ 
dreev reflected as a hole with nearly equal momentum 


[close to kpi for the —>■ process and to kp 2 for 

the —>■ process, as labeled in Fig. 0b)]. Due to 
the momentum matching between the incident electron 
and reflected hole, only zero-momentum Gooper pairs are 
formed within R and the induced pairing amplitude does 
not oscillate as a function of the center-of-mass coordi¬ 
nate of the Gooper pair. The spatial dependence of the 
pairing amplitude can also be inferred from the preser¬ 
vation of time-reversal symmetry; because = E-^^, 
paired states of opposite spin have zero total momentum. 
As shown in Sec. 00 the presence of a triplet pairing 
component again follows from the lifted spin degeneracy. 

Unlike in the ferromagnetic case, the proximity effect 
in 2D S/R junctions is qualitatively different than in ID. 
In 2D, the Rashba spin-orbit vector lies in the plane: 
g = {iky,—ikx,0). Spin-orbit coupling again splits the 
Fermi surface into two Rashba subbands, see Fig. 0c), 
but the crucial difference compared to ID is that these 
subbands can no longer be associated with a definite spin. 
Therefore, for the case where Ep > 0 and \ky\ < kpi, as 
illustrated in Fig. 0c), interband reflections (ei —>■ ^2 or 
€2 ^ hi) are allowed. Because the momentum of an An¬ 
dreev reflected hole within the opposite band differs from 
that of the incident electron, the Gooper pairs that are 
formed by processes of this type have a finite center-of- 
mass momentum. The induced pairing amplitude in two 
dimensions is thus a sum of oscillatory (from interband 
processes) and nonoscillatory (from intraband processes) 
terms. The proximity effect is qualitatively similar when 
Ep < 0; even though only a single Rashba subband is 
occupied, the inner and outer radii of the annular Fermi 
surface play the same role as the two distinct Fermi mo¬ 
menta in Fig. 0c). 


III. SOLUTIONS OF GOR’KOV EQUATIONS 


A. Model and Equations 


We will now show how the physics described in Sec. 01 
follows from the microscopic theory. We consider a two- 
dimensional model of a S/R proximity junction (Fig. 0 
and allow the mass m{x), the SOG constant A(x), the 
Fermi energy Ep{x), and the pairing potential A(x) 
to vary in a stepwise manner across the S/R inter¬ 
face. Specifically, we take m{x) = mpi6{x) + ms0{—x), 
Ep{x) = Epp6{x) + Eps9{—x), \{x) = X0{x), and 
A(x) = A6{—x)i(7y. The model can be represented by 
an explicitly Hermitian Hamiltonian, 


H= 


tp\x) 


-'2^^ 


l{x) 


k^ 

" - Ep{x) 


2m(x) 

tpix) 


^o'y{^X{x)dx A f^x[A(x)]j’ d xX{x)ky 

+ \ (^'ilj\x)A{x)[ilj\x)f A [V'(a::)]^A'l'(x)V'(a;)) |. 


( 6 ) 
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Interfacial scattering is incorporated through a mismatch 
in Fermi velocities/momenta across the S/R interface. 
The Gor’kov equations in this model are given bjl^ 


iOJ + - 9; 


1 „ / 1 


2'"“' V / 2TO(a;) (7x^{x)ky 


dx - 


-lUJ +1^0., 


\ P i 

~ 2m{x) ^ ~ + i9£c[A(x)]} + ax\{x)ky 


Goj,ky (x, x') + A{x)P,y^ky (a:, x') =5{x-x'), 

(7a) 

- A\x)g,y^ky{x,x') = 0. 

(7b) 


We solve the fully quantum-mechanical Gor’kov equa¬ 
tions rather than the quasiclassical Eilenberger equations 
so that we can treat the limit of strong SOC, where the 
splitting of Rashba subbands must be taken into account. 

Even without performing a detailed calculation, it is 
evident from Eq. ([^) that SOC generates triplet pairing 
in the proximity region. If we parameterize the pairing 
amplitude by 

= (s + d • a)icTy, (8) 

then the four coupled Gor’kov equations describing the 
pairing induced in R can be written out explicitly as 

(cfl -I- 5^) s -I- 2mRX{kydx - id^dy) = 0, 

{qr + dl) dx + 2mRX{kyS + dxdz) = 0, 

{a*^ + dl) dy - 2mRiX{kyd:^ - dxs) = 0 , 

(ofl -I- dl) dz + 2mRX{ikydy - dxdx) = 0, 

where we define = 2771 ^( 5 )[iw -|- EpRf^g)] — k^. In 

2D {ky ^ 0), all three triplet components are coupled to 
the singlet component through SOC, whereas only the 
component dy is coupled to s in ID (ky = 0). 


B.l. Solutions in superconductor 

We begin by solving the Gor’kov equations within S, 
which can be expressed as 

{al + dl)g^,ky{x,x') + 2msAidyPl^j^^{x,x') = 0 , 

(10a) 

(as + dl)Pl%^(x,x') + 2msAiayg^j,^(x,x') = 0. 

(10b) 

We denote Green’s functions in the superconductor by 

g5,ky{x,x') = guj,ky(x < 0,a;'), (11a) 

< 0,a:'). (lib) 



B. Calculation Details 


We now proceed with detailed solutions of Eq. Q in 
both the superconductor and Rashba metal. Because 
we are ultimately interested in the pairing amplitude in¬ 
duced in R, we fix x' > 0. In the context of Eq. Q, 
this choice for x' simply ensures that the delta-function 
source term on the right-hand side appears only when 
solving for the Green’s functions in R. 



FIG. 2. Geometry of the S/R proximity effect, together with 
the band structure of each material. System is taken to be 
infinite in the y direction in 2D, while ID corresponds to the 
single-channel limit ky = 0. 
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Solving Eq. (10b) for Q, we obtain 


[x.x = l(J., 


+ dl 


V 1 ^ I (.x,x'). (12) 

2msA ' ’ ” 


Substituting this expression for Q into Eq. (10a), we ob¬ 
tain a fourth-order equation describing 

[(4 + dl){af + dl) + AmlA^]Pl^j^^{x,x') = 0. (13) 


Solving Eq. (|13|) and keeping only those solutions that 
decay in the 


imit x —)■ — 00 , we obtain 


= 5i(^')e-*^" + 55(x')e*^*", (14) 


The Green’s functions can be uniquely expressed as the 
sum of a particular solution to Eq. (19) and the gen¬ 


eral solution to the corresponding homogeneous system 
of equations. 


The particular solution to Eq. (19) is equal to the 


Green’s function of a bulk Rashba metal, which we de¬ 
note by — x'), and is a function of only the dif¬ 

ference X — x'. This function can be obtained by Fourier 
transforming to momentum space, 

(a| - kl)gf^^{K) - 2mR\{ky + ik,)Qf^^^{k,) = 1, 

(20a) 

(4 - kl)gf^^{k,) - 2mR\{ky - ik,)gf^^^{k,) = 0. 

(20b) 


where = 2ms{i^ + EFs) — and 12^ = -I-A^. Here 

and throughout the rest of the calculation, we choose the 
branch cut of the square root function to lie along the 
negative real axis. Note that each coefficient Ci is a 2 x 2 
matrix in spin space and thus contains four unknown 
quantities. 


Solving the algebraic system (20), we find 

2TOit(4 “ML 


-,(0)R 


{kx) = 


g‘'^^{kx) = 


(4 - kl)^ - 4m^A2(A:2 -p k^) ’ 
Amj^X{ky — ikx) 

(4 - kl)'^ - Am\X^{kl + kl)' 


(21a) 

(21b) 


Cl (a;') 


Cl (a;') C 2 (a;') \ 
c?,{x') C 4 (a:') j ■ 


(15) 


Substituting solution ( [l4| into Eq. ( [l^ gives the normal 
Green’s function. 


g^ 


a;,/Ci 


(a:, x') = —it 


Ci{x')e-^P^E^ 


C5{x')t 


where we define rj = cos ^{iuj/A). 


(16) 


B.2. Solutions in Rashba metal 

The Gor’kov equation describing the normal Green’s 
function in R is given by 

[aR + dl - 2mRX{iaydx + axky)]g^,, (x, x') = S{x-x'). 

(17) 

(Note that g^ denotes the Matsubara Green’s function 
in the Rashba metal, 

GS,kyix,x') = g^j^kyix > 0,x'), (18) 

and should not be confused for a retarded Green’s func¬ 
tion.) Equation 0 consists of four equations describing 
the four spin components of the matrix g. This system 
can be easily solved by considering two equations at a 
time; for example (we suppress explicit reference to the 
dependence of the Green’s function on ui and ky here), 


To obtain the Green’s function in coordinate space, we 
perform an inverse Fourier transform. 


The integrals evaluate to 


g^^'>^(x — x') — — 

t ,iD I ulD 

“r /v 2 


iruR 


\x — x' I 


kl 


( 22 ) 


(23a) 


•^2^2 ^ik2S2\x—x'\ 


■^{0)R 


(a; — x') = 


k2 

irriR 

htD I ulD 


X e 


ikisi \x—x' I 


i sgn(a: — a;') — 
i sgn(a; — x') — 


siky 


kl 

S2ky 


X e 


ik2S2\x — x'\ 


(23b) 


Whereas the states in the superconductor are character¬ 
ized by a single momentum p, in the Rashba metal we 
must define two different momenta due to the band split¬ 
ting by SOC, 

ki{2.) = [^1(2)] ~ k^, (24) 

where 

k\^2) = \J + 2,mR{iuj + Err) T R^rX (25) 

are the corresponding momenta in ID. We also denote 
Sa = sgn[lm(fcQ,)], which evaluates to 


(a|. -I- dl)g^^{x,x') - 2mRX{ky + dx)gi^{x,x') = S{x - a:'), 

(19a) 

{qr + dl)gll^{x, x ) - 2mRX{ky - dx)gl\{x, x) = 0. 

(19b) 


f sgn(w) if Err > 0, 

Si = -sgn(w) if Err < 0,a;^ < -AEsoErr, 
I sgn(w) if Err > -AEsoErr, 


S2 = Sgn(a;). 


(26a) 


(26b) 
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In Eq. (26a), we define the spin-orbit energy Eso = 
to_rA^/ 2 (see also Fig. [^. Based on the symmetry of 
Eq. (17), the remaining spin components of the bulk 
Rashba solution can be obtained directly from Eq. (23) 


by noting that Q^, and Gu upon flipping 

the signs of both A and ky. The full bulk solution is given 

by 




r i,iD 


j- 


fci 


UlD I ulD 

/i/“F ii/‘2 

ikisi\x-x'\ _ 'i-rnRS2 

h^D I ulD 


-Si sgn(a: - x')ay 

Ki 


r^i-D 

^2 


-I- S 2 sgn(x - x')ay - 


„ik2S2\x-x'\ 


(27) 


given by 
G5,ky(xix') = 


' f^l|^C9(x') Cio(x') 

C9(x0 -^'i|^Cio(x') ^ 

_i^2||^cn(x') Ci2(x') 

/ f\ ik 2 S 2 — ky ( f\ 

Cll{x') .lO ^ Ci2(r) 

^2 

7(0)H / 


^ik2S2X 


. 

Turning now to the solution for the pairing amplitude 
in the Rashba metal, the Gor’kov equations that we must 
solve are given by 

[ofl +dl- 2mR\{iayda, - a^ky)] (x, x ) = 0. (31) 


Noting that Eq. (31) becomes Eq. 0 upon flipping the 
signs of both uj and ky (save for the delta-function term), 
we obtain the pairing amplitude directly from Eq. (30) 
by making these changes. 


We now seek solutions to the homogeneous system of 
equations 


(a|j -I- dl)G?Ax, x') - 2mR\{ky + dx)Gpx, x') = 0 


{rr + dl)Gi^{x, x) - 2mRX{ky - dj;)GP^{x, x) = 0. 


(28a) 
: 0 . 
(28b) 


Because this is a linear system, it can easily be solved by 
matrix methods. We transform the system (28) into a 
first-order matrix equation. 


•^11, = 


_ »fci^.si + fc„ Ci3(x') Ci4(x') 

ik^si—k. 


Cisix') 


fcfo 

P?iP'‘ cPx') Cieix') 


Ci4(x') 


Ci5(x') 


'lk-2 S2 ky 
]^*1D 


ci6(a;') 


^—ik^six 


^ — ik2S2X 


(32) 


B.3. Enforcing boundary conditions 


Now that we have obtained general solutions to 
the Gor’kov equations in both the superconductor and 
Rashba metal, the sixteen unknown coefficients must be 
determined by boundary conditions. The boundary con¬ 
ditions can be obtained by direct integration of Eq. 0 
over a narrow region near x = 0; they are 


(33a) 

(33b) 




5.x = 


0 

0 

—2mR\ky 


0 

1 

0 \ 


0 

0 

1 


2mR\ky 

0 

1 

to 

X, 

—a^R 2mR\ 

0 / 





(29) 

— \G^ 

d 


The 


four eigenvalues of the matrix in Eq. (29) are Eiki and 


±ik 2 , but we choose only the solutions that decay in the 
limit X —> oo. The prefactors of the exponentials in the 
solution are determined by the eigenvectors. Thus, a full 
solution to the Gor’kov equations in the Rashba metal is 


with the same boundary conditions applying for P as 
well. 


Because the bulk solution (27) contains two linearly 


independent terms, we can separate the dependence of 
the coefficients on x' by writing 


Ci(x') = -I- Cia<i 


ik2S2x' 


(34) 


Each of the sixteen boundary conditions then separates 
into two linearly independent parts, thus giving a total of 
32 boundary conditions that must be solved. The pairing 
amplitude in Eq. (32) becomes 
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xk-\ s 1 ~\~ku 
- \*1D Ci34 

Cl3,l 

pHTD Cl5 1 
Cl5,l “ 


Cl4,l 

ikTsi — k-n 

kf^ ^14,1 
Cl6,l 

xkn S2 ky 
- \.1D Ci6,l ^ 




— z/Cn ^23^ 


Due to the matched momenta in the exponentials of the 


first and last terms of Eq. (351, these terms correspond 


to Cooper pairs with zero net momentum; these terms 
do not oscillate as a function of the center-of-mass coor¬ 
dinate of the pair (x + x')l2. Conversely, the terms with 
mismatched momenta correspond to Cooper pairs with a 
finite momentum; these terms are oscillatory. 

To compactify our notation, we can express the singlet 
and triplet parts of the induced pairing amplitude, as 
defined in Eq. ®, as 


(s,d) = 


2 

E 


[/°/3(‘^>^y)4a/3(w,fcy)]e- 




(36) 

The sums in Eq. (36) run over the two Rashba subbands 
and (/°, f) describes a four-vector of 2 x 2 matrices that 
can be directly related to the spatial dependence of the 
induced pairing amplitude. The four elements of the 
newly defined /* matrices correspond to the four terms 
of Eq. (35); for example, the upper diagonal element of 


y is given by the singlet component of the matrix in the 


first term of Eq. (351, ffi = (ci 4 p — ci3q)/2. For this 


reason, we associate the diagonal elements of the /* ma¬ 
trices with nonoscillatory terms of the pairing amplitude 
and off-diagonal elements with oscillatory terms. 


C. Pairing Amplitude in ID 


Solving the boundary conditions in ID, we obtain = 
/^ = 0, while (X and /° oc CTq. Consistent with our 
previous discussion of Andreev reflection, we find that 
the pairing amplitude is a mixture of singlet and triplet 
components and does not oscillate with {x + x')/2. In 
order to simplify the analytic result, we assume that the 
Fermi energy and the spin-orbit energy are the largest 
energy scales, so that oj,A Eso, Epys)- This allows 
us to expand 

p = PF + 'i^/vF, (37a) 

fci = {kFi + iu}/vR)sgn{EFR), (37b) 

= kF 2 + iojjvR, (37c) 

where vf = \/‘2EFs/n^S = PF/'ms and vr = 

\/X'^ + 2EFR/mR are the Fermi velocities of the super¬ 
conductor and Rashba metal, respectively. We also note 
that, in this limit, si = sgn{uj)sgn{EFR)- The sole 


ikisix' 


ikisix' 


ikT si-\-ky 

-—Ci3^2 Ci4^2 

ikT si — ky 

Cl3,2 —Ci4^2 

ikoS2-\-ky 

Cl5,2 Ci6,2 

^2 

zkfy S2 ky 

Cl5,2 -—Ci6^2 


^—ik^six ^ik^s^x' 


^—ik'^S2X ^ik2S2x' 


(35) 


nonzero triplet component can then be expressed as 


R y) =_ VFAsgn(uj) _ 

^ 2u^u^|a;| + {vj, + 4)12 sgn(F;^^) 

X _ ^-ikF 2 Sgn(ui){x-x')'^ 

(38) 

As is customary, we can also express the nonlocal pair¬ 
ing amplitude in terms of the center-of-mass coordinate 
and the momentum of the relative motion. In this mixed 
representation (x now denotes the center-of-mass coor¬ 
dinate), the triplet pairing amplitude when Err > 0 is 
given by 


dy{x, k) = 


X 


2-KVFVRAsgn{k) 
2vfVr\uj\ E (4 + 4 )^ 

[<5(ei)-^(6)]0(-M, 


(39) 


where ^ 1 ( 2 ) = VR{\k\ — fc^i( 2 )); a similar expression can 
be written if Err < 0. In both cases, the triplet compo¬ 


nent in Eq. (39) consists of terms localized to each of the 


two split Fermi surfaces. However, this triplet component 
vanishes if integrated over the momentum, as the inte¬ 
grated triplet amplitude is proportional to the difference 
in the densities of states on the two Rashba subbands, 


/ dk 

^dy{x, k) oc Ni{Efr) - N 2 {Efr). 


(40) 


Since the subband densities of states in ID are equal re¬ 
gardless of the position of the Fermi level, Ni{Efr) = 
N 2 {Efr) = I/ttvr, dy{x) vanishes. This result is consis¬ 
tent with that of Ref. na which found no induced local 
(equivalently, momentum-integrated) triplet pairing am¬ 
plitude in ID. [Additional Rashba coupling arising from 
a lateral confining potential in quasi-ID wires makes the 
velocities of spin-split subbands, and thus the densities of 
states, different,!^ but this effect is expected to be small 
and ignored here, as well as in Ref. [T71 ] 


D. Pairing Amplitude in 2D 


In 2D, all four components of the pairing amplitude are 
nonzero. While the explicit analytic expressions are too 
cumbersome to be included here, we find, again taking 
w, A ^ Esoj Erri^s)) that the coefficients of Eq. (34) are 
interrelated as follows: 


Cl4,l — —Cl3,l, Ci6,2 — —Cl5,2j 
Cl4,2 = ~C]^3 2 = ~Ci5p = 


Cl4,lj Ci6,2 S ffi, 
Cl4,2 £ C. 


( 41 ) 























FIG. 3. (Color online) Spatial dependence of proximity-induced (a) singlet and (b) triplet pairing amplitudes, as defined 
in Eq. ( |43| l, in a 2D Rashba proximity junction for various strengths of SOC (^s = uf/A is the coherence length of the 
superconductor). Pairing amplitude is plotted in units of its value in the bulk of the superconductor, (c) Quantitative 
comparison between induced singlet and triplet amplitudes for Eso/Efr = 2 [dotted curves from (a) and (b)], which is a 
reasonable choice for BiTel.^ All pairing amplitudes are plotted with tj/A = 0.8, mn/ms = 0.1, Efr/Efs = 0.05, and 
Efs/^ ~ 200. By the Pauli principle, the triplet component has an odd frequency dependence. 


Given these relations, the matrices in Eq. (361 can be 
expressed as 


f° 

r 

fV 

r 


C14.1 R-e[ci4_2] 

R-e[ci4_2] Ci6,2 

— sin0ici4p —Im[sgn(a;)e*®i 014,2] 

-Im[sgn(a;)e“*®2®s"(‘^)ci4,2] sin02Ci6,2 

cos6»ici4,i Re[e*®i"8"(“)ci4,2] 

-Re[e “®®2 ®®"^‘^^Ci 4 , 2 ] - cos 02Ci6,2 

0 jlm[ci4,2] 

—zlm[ci4,2] 0 


sgn(w), 


(42a) 

(42b) 

(42c) 

(42d) 


where we define 0i(2) = sin“^(fcy/fcFi( 2 ))- We see that 
oc CTy, while the remaining components contain both 
diagonal and off-diagonal terms. We therefore conclude 
that intraband Andreev reflection processes contribute 
to the singlet component and to only those triplet com¬ 
ponents that have d || g, while interband processes con¬ 
tribute to all types of pairing. Additionally, each of the 
pairing components has a definite symmetry with respect 
to fcj,; s and dy are even functions of ky, while dx and 
dz are odd. So, while a single trajectory parameterized 
by ky can produce all four pairing components, only the 
components s and dy are nonzero when averaged over all 
possible trajectories. 


To facilitate a quantitative comparison of the induced 
amplitudes of different symmetry in 2D junctions, we 
define a dimensionless analog to the angular-averaged 
quasiclassical Green’s function, whereby we integrate the 


Gor’kov Green’s function over the momentum. 

Note that we explicitly integrate over ky, while the in¬ 
tegration over kx is done implicitly by setting x = x' 
in the nonlocal solution given in Eq. (35). In addition 


to allowing a quantitative comparison between singlet 
and triplet components, integrating over the momentum 
has the added benefit of picking out the odd-frequency 
triplet terms; this follows directly from the Pauli prin¬ 
ciple, as any triplet pairing amplitude is a sum of an 
even-frequency, odd-momentum component and an odd- 
frequency, even-momentum component. 


We calculate the Green’s function defined in Eq. (43) 


numerically without making any approximations. The 
pairing amplitude at a given uj is determined by four pa¬ 
rameters: TTiRilms, Epfi/Eps, Eso/Epp, and Eps/A. 
Figures [^a) and [^b) illustrate the effects of SOC on 
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the induced singlet and triplet amplitudes, respectively. 
While the magnitude of the pairing amplitude is largely 
determined by the Andreev reflection coefficient, which is 
controlled by all four parameters, we find that Eso/Epji 
is the only parameter that controls the singlet-to-triplet 
ratio. Figure 1^ c) compares the singlet and triplet pair¬ 
ing amplitudes using parameters appropriate for a giant 
Rashba semiconductor BiTel, where Ego/Epp « 2.E3 As 
is seen from the plot, the singlet and triplet amplitudes 
can be comparable in magnitude in a real physical sys¬ 
tem. 

As discussed previously in Sec. |IIB[ the generation of a 
triplet pairing component does not require occupation of 
both Rashba subbands. Figure]^ shows that the induced 
triplet component is qualitatively similar when only a 
single subband is occupied. This observation is also rele¬ 
vant for BiTel, as, in samples studied in Ref. [3HI only the 
lowest Rashba subband is occupied at ambient pressure. 


E. Recovering the Quasiclassical Limit 





FIG. 4. (Color online) Quantitative comparison between 
induced triplet pairing amplitude in two-subband (solid line) 
and one-subband (dotted line) cases. Chosen parameters are 
the same as those in Fig.[^ with only a change in the sign of 
Epu/Eps for the one-subband case. 


We will now show that the induced triplet pairing am¬ 
plitude vanishes to leading order in the quasiclassical 
limit, where SOC is taken to be weak (A <C vp) and 
the splitting of the Fermi surfaces is neglected. If there 
is no Fermi surface mismatch between superconductor 
and Rashba metal {mu = ms and Epp = Eps), the 
system can be described by a single Fermi momentum 
kp and a single angle 0 = sm~^{ky/kp). Expanding the 
coefficients from Eq. (41) to linear order in X/vp, we find 


A 


Cl4,l = Ci6,2 = 


2vp cos0(Q 
Cl4,2 = 0{X'^/Vp), 




Because there is only one Fermi momentum, we can 
rewrite Eq. (36) as 


(s,d) = e 


— ik*sx ^iksx 


^ [flf3{ky,i^)Jap{ky,Uj)], 


(45) 

where k'^ = 2m{iuj + Ep) — ky and s = sgn(tLi). We imme¬ 
diately see from Eqs. (421 that all three triplet compo¬ 
nents vanish to first order in X/vp. The nonzero triplet 
pairing amplitude occurs to order X'^/vp (see also Refs. 
[20land [2T|l and is very small in materials with weak SOC. 
Therefore, in order to achieve a sizable singlet-triplet 
mixing, one needs to study the S/R proximity effect be¬ 
yond the quasiclassical limit. 


IV. CONCLUSIONS 

We provided a qualitative physical argument, based 
on the consideration of possible Andreev reflection pro¬ 
cesses, for the presence of triplet Cooper pairing in a 


Rashba material placed in proximity to a conventional su¬ 
perconductor. We also proved the existence of this triplet 
state by solving the exact Gor’kov equations. Even 
though the triplet state vanishes in ID when integrated 
over the momentum, it can be comparable in magnitude 
to the singlet component in 2D when Ego ~ Epp. 

Because triplet pairing occurs only to order X^jvp in 
the quasiclassical limit, where A is the spin-orbit cou¬ 
pling constant and vp is the Fermi velocity, this work 
is most relevant to materials with large A. For exam¬ 
ple, the spin-orbit constant in the bulk Rashba semicon¬ 
ductor BiTel has been reported as A ~ 1 eV-A, corre¬ 
sponding to a spin-orbit energy of Ego = m*X^/2h? ^ 
100 meV (m* ^ 0.lTOe).l^ More conventional semi¬ 
conductors provide a wide range of spin-orbit coupling 
strengths. While GaAs/AlGaAs quantum wells are 
weak spin-orbit materials with A ^ 1 meV-A,®l other 
semiconductor heterostructures, such as InAs/InAlA^Sl 
and InSb/InAlStP^HSl quantum wells, have been re¬ 
ported to have A ~ 0.1 eV-A. Additionally, Ge/Si 
core/shell,l^ InAs,!^ and InStPSl nanowires can reach 
spin-orbit strengths of A ~ 0.1 — 1 eV-A. With so many 
materials belonging to the strong spin-orbit coupling 
regime, it is imperative to understand the proximity ef¬ 
fect in such materials beyond the quasiclassical approxi¬ 
mation. 

Even though we do not treat disorder here explicitly, 
we can comment on its effect. In the presence of impurity 
scattering, the triplet Cooper pairs in a spin-orbit system 
are subjected to spin relaxation. Therefore, in the dif¬ 
fusive limit, where £ and i tgo {£ is the mean 

free path and igo ~ h?/m*X is the spin-orbit length), the 
deca y lengt h of the odd-frequency triplet component is 
~ \JD!Tgo, where Tgo is the spin relaxation time and D 
is the diffusion coefficient. In the diffusive limit of the 
Dyakonov-Perel mechanism,!^ the spin relaxation time 
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is given by r^o ~ so the wave function of the 

triplet Cooper pairs is expected to decay on the scale set 
by iso- However, in materials with strong spin-orbit cou¬ 
pling for which the induced triplet pairing is relevant, the 
spin-orbit length is comparable to the Fermi wavelength; 
therefore, as long as these materials remain “good met¬ 
als” S> 1), the results of this work should apply. 

Recent experiments have measured an electron mobility 
exceeding 200 000 cm^/V-s in an InSb/InAlSb quantum 
well,l2^ thus demonstrating an ability to fabricate ultra¬ 
clean materials with strong spin-orbit coupling. 
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Appendix A: Triplet pairing amplitude as a consequence of lifted spin degeneracy 


In this appendix, we prove that a nonzero triplet pairing amplitude is induced in any ID proximity junction with 
lifted spin degeneracy. For our purposes here, we assume that spin remains an eigenstate of the Hamiltonian, as is 
the case for both a ferromagnet and a ID Rashba material. By calculating the BdG wave functions describing the 
Andreev reflection of an electron to a hole, we will be able to determine the pairing amplitude and show that the 
triplet component is always nonzero when the spin degeneracy is lifted. 

In any nonsuperconducting material where spin is an eigenstate of the Hamiltonian, the scattering wave function 
describing an electron incident on the superconducting interface can be expressed as 


V'TV (a:, e) = e) + a^{£) 








+ r^{e) 



(Al) 


R( T 'l 

where e) is the wave function of the incident electron, is the momentum of a right-moving (left-moving) 

spin-up (spin-down) electron (hole), a^(^)(e) are the Andreev reflection amplitudes, and r^(^)(e) are the normal 
reflection amplitudes. On the superconducting side, the wave function is given by 


ipsix,s) = ti{e) 





hie) 





hie) 



^-iq-x ^ 


tiie) 



^-iq-x 


(A2) 


where = 2msiEps ± — e^), Uq = il + iVA^ — e‘^je)j2 and Uq = (I — j-s/A^ — e^le)j2 are the usual BCS 

coherence factors, and the Uie) are transmission amplitudes. Boundary conditions to be imposed on the wave functions 


are 


V'Af(o) = i’sio), 

VN'fpNiO) = Vstpsi^)- 


(A3a) 

(A3b) 


The velocity operator is a 4 x 4 matrix given by v = dH/di—idx), where H is the Hamiltonian of the BdG equation 
[Hipix.e) = eipix.e)]. 


We first consider the scattering of an incident spin-up electron. In that case, [^iix, e)]^ = (1,0,0,0)e*^f=“. Imposing 


the boundary conditions gives us a set of four equations (a^ = = 0 for this choice of incident wave function); 

l + r^ = huo + tsvo, “ G) = vpsitiuo - t^vo)-, 

ai = hvo + tsuo, vpiai = vpsihvo - huo) 


(A4) 


where is the Fermi velocity of the spin-up (spin-d own ) band of the nonsuperconducting material and vps is the 

Fermi velocity of the superconductor. In obtaining Eq. (A4|, we expanded in the limit e, A «C Ep'jv(S)- We also must 


consider the scattering of an incident spin-down electron, with wave function given by ['0i(a;,e)]™ = (0, l,0,0)e 
The boundary conditions for this second case are 

l + ri = t 2 Uo - t4Vo, 

= -t2Vo+t4Uo, 


ik^ 


i^Fi(l - r^) = vpsihuo + t4Vo), 
— Vp^a^ = Vpsit2Vo + t4Uo). 


(AS) 
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By solving Eqs. (A4) and (A5), we find that the Andreev and normal reflection amplitudes are related by 


r^(e) = 


ai(e) = - —«!(£): 

VFi 

(Mq - vDivFt'VFl - vjps) + VFsjvFt " ^Fj) _ 
{ul - V^){vFtVFl - vj,g) - VpsivFt “ ^Fl) ‘ 




(A6a) 

(A6b) 


The relations of Eq. (A6) will allow us to prove the existence of a triplet pairing component. 


The pairing amplitude is given in terms of the BdG wave functions Ua{x,e) (electrons) and Va{x,e) (holes) by 

■^1/3(2;,a;',£) - Va{x,e)u*fj{x',e). (A7) 

Therefore, the triplet component of the pairing amplitude with zero spin projection is 


J-}{x,x',e)= - J^^^{x,x',e) + J^j^^{x,x',£) ^ V^{x,£)u'^{x',£) + Vl{x,£)u^{x\£). 


(A8) 


These two terms correspond to the two scattering processes considered above. Substituting the wave functions from 
Eq. (All, we find that 


J^l{x,x',£) ~ a^(e) 














.-iki 




(A9) 


Given the relations in Eq. (A 61 , we can reexpress Eq. (A9) as 


VFl 


j^j(x,x',£) ^ a^(£) I e'-’^Th^e -‘"'t. 


+ at(£)r*(£) 


Lx-iktx' XF^ Sgn(£ - A)(ug - ug)(uFtUF; - vlg) + VFsjvF^ " t^F^) ^-ikt 


XFi Sgn(£ - A)(u5 - V^){vFtXFi - xjpg) - Vfs{xf^ - XFi) 




(AlO) 

Note that the triplet amplitude only vanishes when the two spin bands are degenerate, and vf^ = vpi- 

Therefore, triplet pairing is induced by the proximity effect any time the spin degeneracy is lifted. 
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